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The background and development of an exponential model of atmospheric radio re- 
fractivity, the "C.R.P.L. Exponential Reference Atmosphere" is outlined. A set of ray 
tracings for the model is presented in the form of tables of refraction variables for the com- 
plete range of observed values of surface refractive index. A detailed analysis of the ac- 
curacy of the ray tracing and tabulation methods is made for these tables. The variables 
are presented as numbers between one and ten multiplied by the appropriate power of ten, 
thus maintaining a maximum number of significant figures. The tables may be used for the 
solution of practical refraction problems involving elevation angle errors, range errors, and 
similar quant il ies. 



1. Introduction 

Some exponential models of atmospheric radio 
refractivity have been introduced recently to replace 
the standard 4/3 earth treatment for applications 
involving long-range transmission paths, forward 
scatter predictions, or radar tracking at high alti- 
tudes [1, 2, 3, 4, 5, G,]. 2 One of these models, devel- 
oped at the Central Radio Propagation Laboratory 
and called the "C.R.P.L. Exponential Reference 
Atmosphere" has been adopted for use by the Na- 
tional Bureau of Standards in predictions of refrac- 
tion phenomena. 

In the Monograph the background of the ex- 
ponential reference atmosphere is outlined and a set 
of radio ray tracings presented. The ray tracings 
may be used for the solution of many problems in- 
volving refraction. 

2. Background 

In dealingwith refraction problems it is convenient 
to use the radio refractivity denned by 

N=(n-l)Xl0 6 

rather than the radio refractive index, n. The basic 
predictor of refraction effects used in the exponential 
reference atmosphere is the surface value of refrac- 
tivity at the transmitting point, N s . The reasons 
for this become clear upon examination of the normal, 
or average, vertical JV-structure of the atmosphere. 

There is an observed high correlation between the 
mean value of AN, the drop in TV from the surface 
to a height of 1 km, and mean value of N s . This is 
expressed by 



-AiV=7.32 exp {0.005577 N,} 



(1) 



i Digest of NBS Monograph (in press). 

2 Figures in brackets indicate the literature references at the end of tins paper. 



which has been derived from 6 yr means for 45 U.S 
weather stations representing all types of climatic 
conditions encountered in the United States. Thus 
it is evident that on the average the decay of N 
with height may be expressed as a function of the 
surface value of N, at least over the lower portion 
of the earth's atmosphere which is the critical portion 
when considering refraction effects. 

It is also found that the average refractivity struc- 
ture of the atmosphere follows a nearly exponential 
decay with height, the value of TV becoming in- 
significantly small at altitudes over 50 km. 

For these reasons the exponential reference at- 
mosphere was defined as that family of exponential 
refractivity profiles where the rate of N decay with 
height is uniquely specified by the value of N s in 
accordance with eq (1). The exponential reference 
atmosphere is t bus represented by the following set of 
equations: 

N=N S exp { — c e h] 

c e =ln{NJ(N s +AN)} 

with AJV determined from (1) above. 

It is found that several important advantages 
accrue from this particular choice of model, among 
which are: 

1. The TV-profiles involved prove to be a fairly 
good representation of average ^-structure over the 
first 3 km above the surface, a critical range for re- 
fraction effects at low angles to the horizon. 

2. The resulting relationship between the N- 
gradient near the surface and N s is found to be very 
close to the observed average relationship being 
within ±3 percent over the normally encountered 
range of N s (from 300 to 365 for sea level, temperate 
zone stations) . 

3. The resulting profiles and their associated 
gradients are continuous functions of height and 
are in reasonably good agreement with the actual 
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observed refractivity at all heights below the 
ionosphere. 

An illustration of the effectiveness of the expo- 
nential reference atmosphere in predicting average 
refraction effects is given by figure 1, where values 
of angular ray bending are plotted against height 
for rays starting at zero elevation angle. The agree- 
ment between the average bending for 22 climatically 
diverse conditions and the average of the values pre- 
dicted for each climatic condition by the exponential 
reference atmosphere is seen to be quite good, par- 
ticularly as regards the form of the r versus height 
curve. Note, however, that the bending versus 
height predicted by the 4/3 earth model is markedly 
in disagreement with the observed values for most 
heights. 




r at 70 km 
Mean profiles: 15.12 mr 
Mean exp. ref.-. 15.23 mr 
4/3 earth: 30.25 mr 



6 8 10 12 14 

h-h s IN KILOMETERS 



Figure 1. Comparisons of long-term mean bending with 
(a) the average of the values predicted by the appropriate 
exponential reference atmospheres and with (b) those predicted 
for the 4/3 earth's atmosphere. 

The preceding material is intended only as a back- 
ground for the development of the exponential ref- 
erence atmosphere; a more detailed treatment may 
be found in reference [5]. 

The remainder of the Monograph is devoted to a 
detailed analysis of the ray-tracing theory used to 
calculate refraction variables for the exponential 
reference atmosphere, together with an analysis of 
the actual method of calculation. The extensive 
tables of refraction variables are presented in an 
appendix to the Monograph. 

3. Tables of the Exponential Reference 
Atmosphere 

The tables of refraction variables for the "Expo- 
nential Reference Atmosphere" are preseuted in the 
form of ray tracings for profiles with N s =200.0, 
252.9, 289.6, 313.0, 344.5, 377.2, 404.9, and 450.0. 
The uneven values of N s (252.9, 344.5, 377.2, 404.9) 
were chosen to have even values of AN (30.0, 50.0, 
60.0, 70.0). The ray tracings are presented for 
initial elevation angles of 0, 0.5, 1, 2, 4, 8, 15, 30, 65, 
100, 200, 400, and 900 milliradians, the refraction 
variables being tabulated at constant heights (of 
the ray path above the surface) consisting of 0.01, 
0.02, 6.05, 0.1, 0.2, 0.305 (1,000 ft), 0.5, 0.7, 1.0, 
1.524 (3,000 ft), 2.0, 3.048 (10,000 ft), 5.0, 7.0, 10.0, 



20.0, 30.48 (100,000 ft), 50.0, 70.0, 90.0, 110.0, 225.0, 
350.0, and 475.0 km. Thus there are 104 ray trac- 
ings, each tabulated at 24 heights. 

The following refraction variables are tabulated 
at each height: the local elevation angle, 6; the dis- 
tance, d, subtended on the earth's surface by the 
ray path from a height of zero to the tabulated height; 
the angular ray bending, r; the elevation angle error, 
e; the slant range, R , (straight line distance from 
origin to the point on the ray path being considered) ; 
the range error, AR, (the difference between the geo- 
metric length of the curved ray path, R=fdB, and 
the slant range); the radar range error, AR e , (the 
difference between the electrical ray path, R e = 
fndR, and the slant range); Ah, the difference 
between the tabulated height and the height of a 
ray calculated with the 4/3 earth assumption at the 
same distance along the earth's surface; and Ad, the 
difference between 6 for the exponential reference 
atmosphere and 6 for the 4/3 earth ray at the same 
height. All entries are given in a form of floating 
point notation where all numbers are between 1 and 
10 multiplied by the appropriate power of 10, and 
the number ±N.NX10 ±m is printed as ±N.N±m, 
where the second sign is the sign of the exponent 
of ten. 

Four refraction quantities which are not listed 
can be easily obtained from the tables: the angular 
distance subtended by the ray path at the center of 
the earth, cj>=(d/r ) where r is the assumed radius of 
the earth at the ray path origin; the difference 
between the elevation angle of the ray path and the 
elevation angle of the slant range "path", 5=r— e; 
the geometric ray path length, R=R J rAR; the 
electrical ray path length, R e =R J rAR e . 

The elevation angles given in the tables are subject 
to a maximum error of ±0.00005 mr; the values of 
r and e are subject to a maximum error of ±0.001 
mr; the values of the different ranges and range 
errors are subject to a maximum error of ±5 m, 
in particular the range errors are assumed to have 
a maximum error of ±0.5 percent or ±0.1 m, which- 
ever is larger. 

An example of the refraction tables contained in 
the appendix to the Monograph is given in the 
following table for the iV s =313 profile, O =O (313 
is the average value of N s determined from nearly 
two million weather observations taken at 45 different 
weather stations throughout the United States) . 

As an example of the calculation of the four 
unlisted quantities previously mentioned, one obtains 
for the height of 1 km: 

^=^eill4= - 0208800l ' adian 

5 — r-e=5.71671 — 2. 9263 = 2. 7904 mr 

R =i? +Ai?= 133.08495 + 0.00018 = 133.08513 km 

R e =Ro+^Re= 133.08495 + 0.03993 = 133.12488 km. 
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Table 1. Example of tabulated refraction variables for 6 o = and N 8 = 313 



Sul- 
fa cc N" 



313.0 
313.0 
313.0 

313. 
313.0 

313.0 
313. 
313.0 
313.0 
313.0 

313.0 
313.0 
313. 
313.0 
313. 

313.0 
313. 

313.0 
313. 
313.0 

313.0 
313. 
313.0 
313.0 



In it 


Height 


theta 




0.0 


0.010 


.0 


.020 


.0 


.050 


.0 


.100 


.0 


.200 


.0 


.305 


.0 


.500 


.0 


.700 


.0 


1.0C0 


.0 


1.524 


.0 


2. 000 


.0 


3.048 


.0 


5. 000 


.0 


7. 000 


.0 


10. 000 


.0 


20. 000 


.0 


30. 480 


.0 


50. 000 


.0 


70. 000 


.0 


90. 000 


.0 


110. 000 


.0 


250. 000 


.0 


350. 000 


.0 


475. 000 



Delta H 



-5. 1424-4 
-1.0273-3 
-2. 5545-3 
-5. 0557-3 
-9. 9458-3 

-1.4835-2 
-2. 3345-2 
-3. 1325-2 
-4. 2035-2 
-5.6544-2 

-6. 5784-2 

-7. 3319-2 

-4. 6267-2 

2. 4347-2 

1.9795-1 

1. 1676 
2. 5569 
5. 6578 
9. 2212 
1.3025 1 

1.7003 1 
4.7697 1 
7. L613 1 
1.0321 2 



Theta 



1. 49613 
2.11617 
3. 34743 
4. 73735 
6. 70910 

8. 29831 
1.06519 1 
1.26373 1 
1.51633 1 
1.88405 1 

2. 17029 1 
2.70927 1 
3.53155 1 
4.23994 1 
5.15471 1 

7.53020 1 
9.43972 1 
1.22337 2 
1.45417 2 
1.65210 2 

1.82787 2 
2.74520 2 
3.23154 2 
3.73844 2 



Delta 
theta 



7. 6990-2 
1. 0849-1 
1. 6960-1 
2. 3521-1 
3. 2C35-1 

3. 7788-1 
4. 4631-1 
4. 8205-1 
4. 9670-1 
4.4550-1 

3. 4521-1 
7. 6253-3 
-8. 6J62-1 
-1.8936 
-3. 5381 

-8. 7861 
-1.3439 1 
-2.0398 1 
-2.6167 1 
-3.1115 1 

-3.5509 1 

-5.8443 1 

-7.0001 1 

-8.3274 1 



Distance 



1.33678 1 

1.89044 1 

2. 98865 1 

4.22550 1 

5.97338 1 

7.37272 1 

9.43092 1 
1. 11483 2 
1.33073 2 
1.63885 2 

1. 87356 2 
2. 30289 2 
2. 92785 2 
3. 44203 2 
4. 07963 2 

5. 65227 2 

6. 87953 2 
8.66238 2 
1.01334 3 
1. 13949 3 

1.25151 3 
1.83614 3 
2. 14610 3 

2. ir,9i5 :; 



6. 01363-1 

8. 50046-1 

1. 34194 
1.89273 

2. 66351 

3. 26995 
4. 14574 
4. 85517 
5.71671 

6. 87418 

7. 69459 
9. 04109 
1.06243 1 
1. 16082 1 
1.24648 1 

1.33857 1 
1.35470 1 
1.35807 I 

1.35824 1 
1. 35825 1 

1.35825 1 
1.35825 1 
L. 35825 ! 
1. 35825 1 



Error 


angle 


3. 0068-1 


4. 2516-1 


6.7171-1 


9. 4849-1 


1. 3382 


1. 6474 


2. 0975 


2. 4678 


2. 9263 


3. 5600 


4. 0267 


4. 8367 


5. 9039 


6. 6859 


7. 5264 


9. 0531 


9. 8441 


1. 0615 1 


1.1051 1 


1. 1335 1 


1. 1540 1 


1.2210 1 


1.2419 1 


1.2584 1 



Slant range 



1.3367819 1 

L. Si 104449 1 

2.9886667 1 

4.2255465 I 

5.9734900 1 

7. 3729190 1 
9.4313420 1 

1.1149053 2 
1.3308495 2 
1.6390785 2 

1. 8739027 2 
2. 3035214 2 
2. 9291690 2 
3. 4442178 2 
4. 0833598 2 



5. 6628165 2 

6. 8993595 2 

8. 7039820 2 

1. 0202252 3 

1. 1495090 3 

1.2650293 3 

1.8820256 3 

2. 2215878 3 

2. :.S7. r .l. r .3 'A 



0.01)0 
1. 000-6 
3.000 6 
6.000 6 

2.000-5 

3. 500 •". 

6. 600-5 
1.00(1 l 
1. 800-4 
3. 100-4 

4. 600-4 

7. 800-4 

1. 400-3 

2. 000-3 
2. 790-3 

4. 460-3 

5. 470-3 

6. 520-3 
7. 100-3 

7. 500-3 

7. 800-3 

8. 690-3 
8. 900-3 
9. 000-3 



Delta 
R /: 



I. IM 3 

5. 912-3 
9. 334-3 
1.317-2 
1. 854-2 

2.27S 2 
2.890-2 

3. 387-2 
3.993 2 

4. 812-2 

5. 398-2 

6. 369-2 

7. 539-2 

8. 284-2 

8. 960-2 

9. 801-2 
1. 003-1 
1. 017-1 
1. 024-1 
1. 028-1 

1. 031-1 
1. 038-1 
1. 042-1 
1. 042-1 
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